Interlayer Magnetic Coupling and the Quantum Hall Effect in Multilayer Electron 

Systems. 
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We study the effect that the electron-electron interaction has on the properties of a multilayer 
electron system. We consider the case corresponding to filling factor unity in each layer. We find 
that as a function of the sample parameters the system has ferromagnetic, canted antiferromagnetic 
or paramagnetic interlayer spin correlations. These three ground states are QHE phases, because of 
the existence of a finite activation energy. In the ferromagnetic phase the gap is due to the intrawell 
exchange energy, whereas in the paramagnetic phase the gap appears due to the spatial modulation 
of the inter well coherence. 



PACS number 73.40.Hm 

The quantum Hall effect (QHE) is one of the most 
striking phenomena observed in two dimensional electron 
gas (2DEG) systems 0], The QHE occurs because the 
2DEG becomes incompressible at certain filling factors. 
In the fractional QHE, the energy gap source of the in- 
compressibility is produced by interactions between elec- 
trons. In the even integer QHE the incompressibility is 
due to the quantization of the electron kinetic energy and 
interactions are not important for the occurrence of the 
QHE. In the case of the odd integer QHE the incompress- 
ibility occurs because at odd filling factors the 2DEG is 
ferromagnetic. Due to the electron-electron interaction 
the ground state of a 2DEG at odd filling , is completely 
spin polarized even for zero Zeeman energy [?]. Given 
the new physics which appears in 2DEG in the QHE 
regime, the question that arises is whether the quantum 
Hall phases are unique to 2DEG or they can occur in 
three dimensional (3D) conductors ||. In this direction, 
some studies in narrow gap 3D semiconductors in the 
strong magnetic field limit have shown some signatures 
of an incipient quantum Hall phase S. On the other 
hand the progress in epitaxial growth has made possible 
to fabricate semiconductor systems where 2DEG's with 
extra degrees of freedom exists: 

i) Wide parabolic quantum wells where a thick electron 
gas layer (~2000A) is formed. This system present a 
clear QHE phase [Q. 

ii) Double layer 2DEG systems with electrons confined to 
two parallel sheets separated by a distance comparable to 
that between electrons within a plane ||. Double layer 
systems present QHE at total integer filling factors, even 
in absence of tunneling between the electron planes. 

iii) Superlattices, where an appreciable dispersion of the 
electronic spectrum in the direction perpendicular to the 
layers exits. Accurately quantized Hall plateaus have 
been observed in these multilayer systems [[U?]] when a 
magnetic field is applied parallel to the superlattice axis. 
Studies of vertical transport in these supperlattices have 
shown Q the existence of a chiral two-dimensional sys- 
tem that form at the surface o the layered system || . 
The poor mobility of the bulk narrow gap semiconduc- 
tors and the small number of new degrees of freedom of 



the parabolic and double quantum wells with respect the 
2DEG, made the superlattices the best candidates for 
studying QHE phases in 3D conductors. 
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FIG. 1. Phase diagram for a multilayer system, with 
filling factor unity in each well. The Zeeman coupling is 
H=0.01e 2 /e£ and the electron layer thickness b — 0.8£. Three 
phases are present: a ferromagnetic phase (shadow region), 
a canted antiferromagnetic region (C) and a paramagnetic 
region (P). 

In this work we study the effect that the electron elec- 
tron interaction has on the properties of the multilayer 
electron system. We consider the case corresponding to 
filling factor unity in each electron layer. Two points are 
raised in this paper: the magnetic order of the electron 
layers and the conditions for the occurrence of the QHE 
in this system. The main results we obtain are the follow- 
ing: i)As a function of the sample parameters (Zeeman 
coupling, H, interlayer tunneling, t, and barrier thick- 
ness, d) we obtain that the system changes from a QHE 
state with interlayer ferromagnetic spin correlations to a 
new QHE state with canted antiferromagnetic interlayer 
correlations (see Fig.l). For larger values of the tunneling 
amplitude, we obtain that the system undergoes another 
phase transition towards a paramagnetic state. These 
transitions are second order phase transitions, ii) We also 
study the value of the activation energy as a function of 
the sample parameters. In Fig. 2 we plot this energy gap 



1 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



t (e 2 el) 

FIG. 2. Activation energy of a multilayer system as a 
function of the tunneling amplitude. The Zeeman coupling is 
H=0.01e 2 /e£, the barrier thickness d = 0.2£ and the electron 
layer thickness is b = 0.8£. The vertical lines indicate the 
values of t where the different phase transitions occur. 

as a function of the tunneling amplitude for a multilayer 
system with d = 0.2£ and H = O.Ole 2 /^ (here I is the 
magnetic length). We find that this gap is finite even for 
very large values of t, where the system is paramagnetic. 
This implies that the paramagnetic phase of the multi- 
layer system is also a QHE phase. As we explain below, 
the energy gap in the paramagnetic phase appears be- 
cause the system breaks spontaneously the translational 
symmetry along the multilayer axis by modulating the 
interwell coherence. From our results we conclude that 
in multilayer electron systems, with filling factor unity 
in each well, the QHE prevails in all phases due to the 
existence of a finite activation energy, even at d, t — ► oo. 

We treat the electron electron interaction in the 
Hartree-Fock (HF) approximation. In 2DEG's in the 
QHE regime, this approximation || gives results which 
agree well with experiments fLOf . In double quantum 
well systems, at total filling factor unity, the HF approx- 
imation is less accurate because of the competition be- 
tween inter and intrawell correlations, nevertheless for 
small values of the layer separation the agreement with 
experiments is rather good || . In double layer systems at 
total filling factor two, the HF approximation gives rea- 
sonable results because of the existence of filled Landau 
levels and the existence of charge excitation gaps, which 
make the HF ground state to be a good approximation to 
the real many-body ground state ]TT[ |. From the above, 
the HF approximation is expected to be a good approxi- 
mation for describing weakly coupled multilayer systems 
at filling factor unity in each well. For stronger cou- 
pled multilayer systems we expect HF to be a reasonable 
approximation for the real many-body ground state, be- 
cause the existence of filled bands and charge excitation 
gaps. 



The calculations presented here employ realistic 
Coulomb interaction potentials, take into account inter- 
layer tunneling and Zeeman coupling and therefore we 
expect our results to be qualitatively and quantitatively 
trustworthy. We take the multilayer vertical axis as the 
^-direction and the electrons live in the x — y planes. The 
magnetic field, B, is applied in the z-direction. Since B is 
very strong, we only consider states in the lowest energy 
Landau level of the lowest energy subband of each well. 
The Hamiltonian of the system is written as H — Ho + V, 
with 

' E { C ta, k C^,k + h.c) 

<i,j> ,k,a 

- H Y. aC t,*Ci,<r,k , (1) 

where Cf a k creates an electron in the lowest Landau 
level in layer i with spin a (a = ±1) in the direction 
of B, and with intra Landau level index k. The sum in 
the first term of Eq.l is over first neighbors layers. The 
many-body part of H takes the form, 

v = ^EE E v^e-^/v^-*')^ 

(7,(7' i,j fc,/c',q 
X ^i,a,k+q y ^j,(j' ,k'Cj,v' ,k'+q y Ci,(7,k , (2) 

where S is the sample area and the interaction potential 
has the form Vi.j — 2ire 2 /eq Fij (q, b, d), with Ffj being 
finite layer thickness form factors | fL2| , which depends on 
the thickness of the 2DEG in each layer, b, and on the 
barrier thickness, d. 

The multilayer systems are doped in the barriers and 
in absence of B, all layers have the same number of elec- 
trons. We consider that this is also the situation at large 
values of B, since any other situation it would cost a large 
Hartree energy. At filling factor unity in each layer the 
interwell correlation is not very important so we only con- 
sider solutions with translational symmetry in the plane 
[x, y) of the electron gases. Broken translational symme- 
tries along the multilayer axis (z-direction) are allowed, 
always with the condition of having filling factor unity 
in each layer. With these constrains, the HF expectation 
value of V takes the form, 

<^>=-^EE^' ¥/2 

< ^VT,fcCj>',fc >< ^t,a' .k-qy^i^M-qy > , (3) 

here the sum in k is over all its possible values. By min- 
imizing the energy < H >=< H > + < V >, we obtain 
the energy of the ground state of the system and its prop- 
erties. 

We have solved the Hamiltonian for different values of 
d, t and H. For each layer i we calculate the expectation 
value of the total spin operator per electron < Si >, 
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Si, x + iSi,y = TT X! ^l+l.fc^*.- 1 .* ' ( 4 ) 

here iV^ = S/2tt£ 2 is the the Landau level degeneracy. 
For characterizing the ground state it is also necessary 
to quantify the interlayer coherence which is given by 
the following expectation value, 

^(i) = ^J2 <C t,kCi + ^,k> ■ (5) 
* k 

This quantity represents the coherence between wells. 
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FIG. 3. Expectation values of the total spin operator per 
electron as a function of the tunneling amplitude. The Zee- 
man coupling is H=G.01e 2 / el, the barrier thickness d — 0.21 
and the electron layer thickness is d = 0.9,1. The vertical tick- 
marks in the lower x axis indicate the values of t where the 
different phase transitions occur. 

Looking to the values of < Si > we find three different 
classes of ground states (see Fig.l and Fig. 3): 

1) Ferromagnetic phase, where all electron layers are fully 
spin polarized in the direction of the magnetic field, i.e. 
< Si >= (0, 0, 1/2). This phase occurs for small values of 
t or large values of H . In the ferromagnetic phase the in- 
tralayer coherence is more important than the interlayer 
coherence and all the expectation values of the operators 
A ff)<r /(i) are zero, and there is not vertical kinetic en- 
ergy contribution to the total energy. The ferromagnetic 
ground state is a QHE phase, and the activation gap (see 
Fig. 2) is the cost in energy of adding an electron to the 
system with the spin pointing antiparallel to the mag- 
netic field. In this phase the ground state has the same 
translational symmetry than the original Hamiltonian. 

2) Canted antiferromagnetic phase. In this phase the total 
spin in each layer acquires a component perpendicular to 
the magnetic field, < Si >= (S^j., Si >g ), and the magni- 
tude of < Si > is smaller than its maximum value, 1/2. 
In this phase the sign of Sj ; j_ alternates from layer to 



layer, i.e. Si.j_ = — S,±i so that the translational sym- 
metry of the Hamiltonian is spontaneally broken and the 
unit cell in the z-direction consists of two electron layers 
which are labelled 2i — 1 and 2i. The z-component of 
< Si > is finite and it has the same value in all layers, 
therefore in the canted phase there is an interlayer an- 
tiferromagnetic coupling of the transverse component of 
the total spin |l3| ] . By performing calculations in bigger 
size unit cells, we have checked that this phase is stable 
with respect spiral ordering of the transverse component 
of < Si >. 
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FIG. 4. Interlayer coherence parameters as a function of 
the tunneling amplitude. A ++ (i) are real and we plot its real 
part. A+_(i) are imaginary and we plot its imaginary part. 
The Zeeman coupling is H=0.0le 2 /el, the barrier thickness 
d — 0.21 and the electron layer thickness is b — 0.81. The 
vertical dashed tickmarks in the upper x axis indicate the 
values of t where the different phase transitions occur. 

In the canted phase the interlayer coherence parameter 
is different from zero and verify the relations: 

A +1+ (i) = A_,_(i) 

A + ,_(i) = -A*,+(i) , (6) 

but the interlayer coherence parameter depends on i. 
In Fig.4 we plot A +j+ (2i), A +!+ (2i - 1), A+,_(2i) and 
A + ,_(2z — 1), as a function of t. We see that 

A CT)CT ,(2i) ? A fflff /(2*-l) , (7) 

and this implies that there is a modulation of the in- 
terlayer coherence. Therefore in the canted phase the 
translational symmetry along the multilayer axis is bro- 
ken not just by the antiferromagnetic ordering of the lay- 
ers, but also by the modulation of the interlayer coher- 
ence. In the selfenergy calculation, the modulation of the 
interlayer coherence acts as a spatial modulation of the 
hopping amplitude, and this modulation contributes to 
the opening of an energy gap at the Fermi level. The 
canted phase appears at intermediate values of the tun- 
neling amplitude, see Fig.l and Fig. 3, and the reason 
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for its existence is that in this phase the system can take 
advantage of the kinetic energy by creating interlayer an- 
tiferromagnetic spin correlations. The antiferromagnetic 
order is canted in order to minimize the lost of Zeeman 
energy. The canted ground state is a QHE phase. The 
transport activation energy is finite because the system is 
partially spin polarized in the direction of B, and because 
the interlayer coherence is spatially modulated. 

Canted ground states corresponding to rotations of all 
the Si y ± are degenerated and therefore this phase should 
get a gapless collective mode associated with this degen- 
eracy. 

3) Paramagnetic phase In this phase the expectation value 
of the total spin operator is zero in all layers, < >= 0. 
This phase occurs at large values of t, where the ki- 
netic energy and interwell coherence energy is much big- 
ger than the Zeeman and intrawell exchange energy. In 
this phase A+ _(i) = but the equal spin interwell co- 
herence parameters arc different from zero and verify 
A+,+ (2i - 1)=A_ _(2i - l)^A+,+ (i)=A_ _(2i). In this 
phase, the system breaks spontaneously the translational 
symmetry by modulating the interwell coherence along 
the vertical axis of the multilayer. In the paramagnetic 
ground state the unit cell consists of two electron lay- 
ers. This modulation of the interwell coherence creates 
an energy gap at the Fermi energy, and the paramagnetic 
ground state is a QHE phase. 

The superlattices studied in references ||,[7| have thick 
barriers and they are in the ferromagnetic phase. For 
studying antiferromagnetic and paramagnetic phases, it 
is necessary multilayers with thin barriers and large tun- 
neling amplitudes. The multilayers are usually doped in 
the barriers and in the case of thin barriers this produces 
a strong scattering of the electrons by impurities, which 
prevents many-body driven ground states. It is possible 
to circumvent this problem by working with superlattices 
superimposed on wide parabolic wells. These 
systems are remotely doped and it is possible to obtain 
multilayers with thin barriers and high electron mobility. 
The ground states magnetic properties can be studied 
experimentally by using optically pumped nuclear mag- 
netic resonance. This technique has been very useful for 
the study of the magnetic nature of 2DEG's jl7). Also it 
could be very useful the application of a magnetic field, 
Bu , parallel to the electron sheets. Bu changes the value 
of the Zeeman coupling, and as a function of its strength 
the ground state of the system could change. This phase 
transition could be identified by studying the activation 
energy as a function of Bu [Q. Strong enough Bu also 
destroys the interlayer coherence In the param- 

agnetic phase the activation energy is due to the spatial 
modulation of the interlayer coherence, and the applica- 
tion of a strong Bu , would destroy the QHE. 

In conclusion, we have studied the effect that the 
electron-electron interaction has on the properties of a 
multilayer electron system. We consider the case corre- 



sponding to filling factor unity in each layer. We have 
found that as a function of the sample parameters the 
system has ferromagnetic, canted antiferromagnetic or 
paramagnetic interlayer spin correlations. We have ob- 
tained that these three ground states are QHE phases, 
because of the existence of a finite activation energy. In 
the ferromagnetic phase the gap is due to the intrawell 
exchange energy, whereas in the paramagnetic phase the 
gap appears due to the spatial modulation of the inter- 
well coherence. 
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